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We investigate the possibility of ultracold neutron (UCN) storage in quantum states defined by 
the combined potentials of the Earth’s gravity and the neutron optical repulsion by a horizontal 
surface of liquid helium. We analyse the stability of the lowest quantum state, which is most 
susceptible to perturbations due to surface excitations, against scattering by helium atoms in the 
vapor and by excitations of the liquid, comprised of ripplons, phonons and surfons. This is an 
unusual scattering problem since the kinetic energy of the neutron parallel to the surface may be 
much greater than the binding energies perpendicular. The total scattering time constant of these 
UCNs at 0.7 K is found to exceed one hour, and rapidly increasing with decreasing temperature. 
Such low scattering rates should enable high-precision measurements of the scheme of discrete energy 
levels, thus providing improved access to short-range gravity. The system might also be useful for 
neutron /3-decay experiments. We also sketch new experimental concepts for level population and 
trapping of UCNs above a fiat horizontal mirror. 


I. INTRODUCTION 


Slow neutrons play an important role in low-energy par¬ 
ticle physics as a tool and an object, in investigations 
of the free neutron’s properties and its interactions with 
known or hypothetic fields with high precision [IHi- A 
particular class of experiments employs neutrons with en¬ 
ergy lower than the neutron optical potential of typical 
materials, i.e. in the order of up to 300 neV. These so- 
called ultracold neutrons (UCNs) can be emprisoned for 
many hundreds of seconds in well-designed ’’neutron bot¬ 
tles” . By virtue of the neutron magnetic moment of 60 
neV/T magnetic trapping is feasible too, and also the 
gravitational interaction with a potential difference of 100 
neV per meter rise can play its role in UCN storage and 
manipulation Si- Perhaps the most prominent appli¬ 
cation of UCNs is the search for a non-vanishing electric 
dipole moment of the neutron. Finding a finite value or 
just improving current limits is a way of investigating new 
mechanisms of CP violation beyond the standard model’s 
complex phase of the weak quark mixing CKM matrix, 
and the matter-antimatter asymmetry in the universe 7|. 
Latest experimental results were published in Refs. [^, 3- 
Similarly long standing are the experimental efforts to de¬ 
termine the neutron lifetime with high accuracy. Its value 
enters the calculations of weak reaction rates in big-bang 
nucleo-synthesis and stellar fusion [13, El , and it is also 
crucial for deriving the weak axial-vector and vector cou¬ 
pling constants of the nucleon needed to calculate many 
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important semi-leptonic cross sections. 


Discrete energy levels of UCNs in the Earth’s gravi¬ 
tational field were proposed by Lushikov and Frank in 
1978 and demonstrated experimentally in the past 
decade liml. Precise measurements of the energy lev¬ 
els, that are not equi-distant, offer an interesting tool for 
tests of various new scenarios of particle physics. The 
range of effects investigated is determined by the charac¬ 
teristic size, several tens of micrometers, with which the 
neutron wave functions are bound in one dimension. De¬ 
viations from the Newton’s gravity law at small distances 
can for instance be interpreted as a signal of large extra 
dimensions at the sub-millimeter sc ale IlfiL 17| or as a hint 
for dark-energy ” chameleon” fields [l3 [l3] • A recent de¬ 
velopment, called gravity resonance spectroscopy (GRS), 
where transitions between levels are induced by vibrating 
the mirror, has paved a way towards sensitive tests of such 
scenarios [20l| (note, however, a strong competition from 
atomic physics in the chameleon search [Ej). A compet¬ 
ing, alternative, method will employ oscillating magnetic 
field gradients [^, . The GRS experiment described in 

Ref. [la] has already set stringent limits on chameleons. It 
has also constrained axion-like particles, improving the re¬ 
sult of an analysis of the non-resonant gravity experiment 
described in Ref. [2^. A method not relying on spatial 
quantum states of the neutron employs spin precession of 
trapped UCNs close to a heavy mirror [^. The sensitiv¬ 
ity of the GRS experiment was almost as good as a first 
search of that latter type [^ and recently much improved 
[^ , still with potential for large further gains in sensitiv¬ 
ity. For the gravity experiment too, a large gain is still 
to be expected, notably once an adaptation of Ramsey’s 
molecular beam technique of separated oscillatory fields 
to GRS is implemented [28|. In addition, a search for a 
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non-zero neutron charge based on the latter technique has 
been proposed 

All current experiments on gravitational quantum 
states of the neutron employ highly polished quartz mir¬ 
rors. These are expensive, limited to sizes of several tens 
of centimeters, and they have to be horizontally levelled by 
some active means. In this respect, using a liquid surface 
as a mirror might initiate a qualitatively new approach. 
On one hand it may furnish a remedy to the aforemen¬ 
tioned limitations. On the other hand, the interactions 
of neutrons prepared in gravity states with excitations or 
structural decorations of the liquid surface could enlarge 
the applications of these states to investigations of the 
surface physics. The present article provides a theoret¬ 
ical investigation of the possibility to store UCNs in the 
lowest gravitational energy states on the liquid helium sur¬ 
face, by analysis of scattering by helium atoms in the gas 
phase and by various excitations in both the bulk and at 
the surface of the liquid helium. Obviously, a long stor¬ 
age time constant is a necessary condition for conducting 
experiments using a mirror made of this quantum liquid. 
A separate section sketches some experimental concepts 
addressing issues arising in real studies employing those 
neutrons, notably population, trapping and detection. 

Properties of neutrons on the liquid helium surface are, 
in several aspects, similar to those of electrons. The two- 
dimensional electron gas on a surface of dielectric media 
has for many decades been a wide subject of research (for 
reviews see, e.g., (sol - f^ l. In contrast to the gravitational 
force in the neutron case, the electrons are attracted to 
the boundary by the electric image forces through which 
they become localized in the direction perpendicular to 
the surface. The surface of superfluid helium has no solid 
defects (like impurities, dislocations, etc.) and offers a 
unique opportunity to create an extremely pure 2D elec¬ 
tron gas. The mobility of electrons in this gas usually 
exceeds more than thousand times that of electrons in 
2D quantum wells in heterostructures. The system thus 
simulates a solid-state 2D quantum well without disorder. 
Many fundamental properties of a 2D electron gas have 
been studied with the help of electrons on the surfaces of 
liquid helium. Various electronic quantum objects can be 
experimentally realized on the liquid helium surface, such 
as quantum dots [s^, ID electron wires [s^, quantum 
rings and others. The electrons on the liquid helium 
surface may also serve for an experimental realization of a 
set of quantum bits with very long decoherence time . 
If neutrons can be made to rest in surface states in suffi¬ 
cient densities we can hope for comparable studies using 
neutrons rather than electrons, and possible new states of 
quantum matter. 


II. NEUTRONS ABOVE A FLAT HELIUM 
SURFACE 


We consider a plane boundary between superfluid "^He 
(situated at vertical coordinate 2 < 0) and its saturated 
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FIG. 1: Schematic representation of the vertical potential and 
the first two states of a neutron above a horizontal mirror of 
liquid helium. 


vapor (z > 0). The interaction of a neutron with a ^He 
atom with nuclear coordinate Ri can be expressed as a 
Fermi pseudo-potential given by 

U, (r) = (r - R,) = (r _ R^), (1) 

m 

where <5^^^ (r) is the 3D Dirac delta function. Substitut¬ 
ing the bound coherent neutron scattering length ape = 
3.26 X 10“^^ cm of a “^He atom and the neutron mass m = 
1.675 X 10“^^ g, one obtains the value U = 1.36 x 10“'*^ erg 
cm^ [1 erg = 10“^ J ~ 6.25 x 10^^ eV]. Within the helium 
bulk the neutron interacts with a ’’forest” of 5-function 
potentials with volume concentration given by the parti¬ 
cle density of ^He atoms, nue ~ 21.8 nm“^ at T < 2.18 
K. As a result of the interference of a plane incident wave 
with the spherical scattered waves from each ^He nucleus, 
neutron propagation in the bulk can be described by a 
constant neutron optical potential given by the spatially 
averaged pseudo-potentials of the helium atoms in a vol¬ 
ume If containing many atoms, i.e. 

Vo = = Unue ( 2 ) 

Ri 2.965 X 10“^° erg = 18.5 neV. 

Above the ^He surface, neglecting interactions with the 
helium vapor discussed further below, the neutron is ex¬ 
posed to the gravity potential, 

V{z) = mgz, (3) 

where g = 981 cm/s^ is the acceleration of free fall (grav¬ 
itational acceleration at Earth’s surface), giving the grav¬ 
ity force for the neutron, mg = 1.643 x 10“^^ erg/cm = 
1.026 neV/cm. 

One can easily solve the one-dimensional Schrodinger 
equation for a neutron in the potential given by Eqs. m 
and ([ni) and shown in Fig. [TJ The corresponding Hamil¬ 
tonian is given by 

h^A 

Ho = —--h mgz + VoO (-z), (4) 

2m 
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where A = is the Laplace operator, V is the gradient 
operator, and 9 (x) is the step function [9 (x) = 1 for a: > 0 
and 9 (x) =0 for cc < 0]. The x, y and z coordinates 
separate in this equation, and the neutron wave function 
is given by a product 

=V'|| (i’ll) V'xW, (5) 

where r|| = {x,y) is the 2D coordinate vector along the 
surface, while r = (x, y, z) stands for the 3D coordinate 
vector. In the x-y plane, the neutron wave function is 
given by a normalized plane wave, 

V'll i^w) = exp (iP||r||/h), (6) 

where S is the He surface area, and P|| is the 2D neutron 
momentum along the surface. For z < 0 we can neglect 
the weak gravitational potential given by Eq. m as com¬ 
pared to the much stronger potential wall given by Eq. 
@. The z-dependent part of the neutron wave function 
in this region is then approximately given by 

ijjj_ (z) = (0) exp (kz) (z <0, E± < Vq) , (7) 

where k = ik± = 2m (Vq — Ej_ )/ h, and E± is the 
neutron kinetic energy along the z-axis. For E± <C Vq, 
Kq = \/2mVQ/h K, 2.4 X 10® cm“®, i.e. the neutron pene¬ 
tration depth into the liquid helium is 

Kg ^ ~ 33 nm. ( 8 ) 

For z > 0, the neutron wave function is given by 

^x{z) = C Ai[{z-Ex/mg)/ zq] (-z > 0), (9) 

where C is a normalization coefficient, Ai {x) is the Airy 
function, and 

zo = {h^= 5.87 /rm (10) 

is a characteristic length scale of the neutron wave func¬ 
tion in low energy states. For z > E±/mg the wave func¬ 
tion given in Eq. decreases exponentially. For Ex < Vq 
the energy spectrum of the z-axis motion of a neutron 
above the liquid-helium surface is quantized. The eigen¬ 
values of energy are given by the boundary condition at 
z = 0 , i.e. 

^'x (0) _ _ ^/2m{Vo - Ex) _ Ai'{-u) 

tpxi^) ^ zoAi(—u)’ 

where the prime indicates the derivative, and 

u = Ex/mgzQ. ( 12 ) 

Equation (HB) for finite Vq can be solved only numeri¬ 
cally. A characteristic scale of separations between low- 

1 /3 

est energy levels is given by mgzQ = {p?"mg^ /2') = 

0.96 X 10-24 erg = 0.6 peV. 


In the limit Vg —>■ oo the energy levels are given by 

En = mgzoan+ 1 , n = 0,1,..., (13) 

where —a„+i are the zeros of the Airy function, 

ai = 2.338, a2 = 4.088, ag = 5.521, a/^ = 6.787. 

(14) 

For n » 1, 

a„ ~ (37rn/2)^^® — (7r2/96n)^^^ . (15) 

Note that this expression provides a highly accurate ap¬ 
proximation even for small n: for oi\ the error is only 
|Aa/a| ~ 0.0018, and further decreases for n > 1 , e.g., 
jAa/aj ~ 0.0005 for a^- 
For finite Vg, Eq. m can be rewritten as 

Ai'{-u) = Ai{-u) (16) 

using the dimensionless constant 

r^g = Eg/mgzg = 3.078 X 10^, (17) 

evaluated on the r.h.s. for Vg of superfluid He. Eqs. dm) 
and CB give the following values of = En/mgzo of 
the discrete energy of a neutron above the He mirror: 

af® = 2.332, = 4.082, = 5.515, af® = 6.781. 

(18) 

These values are indeed very close to the values for 
Vg —>■ oo [compare Eq. (fCT) ]. because ? 7 q 1. The neutron 

wave functions above liquid He are also very close to those 
for Vg —>■ oo, except for a region close to z = 0, where they 
acquires values ^’xn ( 0 ) = Ai (—q;[[^® g) which are small 
but Hnite, since ~ ctn+i- For the n-th energy level 

the normalization coefficient is given by 


pCO 

/ -r \ 

2 

/ dz 

Ai(--aZi) 

1 1 A i ^ He ^ 1 ^ ^ — 2kz 

lo 

V^o J 



The first three normalization coefficients are given by 
Co ~ 1.4261/y^ ~ 59 cm-^/ 2 ^ (j^ ~ 5^.3 
C 2 ~ 47.8 cm-®A, Using the values Ai (—af®) Ri 0.004, 
Ai(-Q;f®) R! -0.0045, and Ai (-af®) r 0.0048, this 
gives tpxo (0) ~ 0.236 (0) r —0.231 cm-^A^ 

ipX 2 (0) ~ 0.23 cm-^/ 2 ^ These values will be used 
further below. 

In the subsequent sections we consider the stability of 
a neutron in a bound surface state against various scat¬ 
tering processes. Note that we deal here with a rather 
unfamiliar scattering problem in that the kinetic energy 
of the neutron parallel to the surface may be many or¬ 
ders of magnitude greater than the binding energies in the 
perpendicular direction. We calculate the temperature- 
dependent scattering rates Wvap, Wrip and Wsur due to 4He 
atoms in the vapour above the surface, due to waves on 
the helium surface, called ripplons, and due to excitations 
called surfons, respectively. 














4 


III. SCATTERING OF NEUTRONS IN 
SURFACE STATES BY HELIUM VAPOR 


^He vapor atoms can be considered as point-like im¬ 
purities with interaction potential given by Eq. ([T]). The 
momentum distribution of the vapor atoms is given by the 
Bose distribution [s^ . 


Np = exp 


f — Eae \ 

V kBT J ’ 


( 20 ) 


where fee = 1-38 x 10“^® erg/K is the Boltzmann con¬ 
stant, /^/ks = —7.17 K is the chemical potential of liquid 
'*He (evaporation energy of a "'^He atom) for T —?> 0, and 
Bue = j2M is the kinetic energy of a "'^He atom with 
momentum P. The total atom concentration in the vapor 
is given by the integrated momentum distribution, i.e. 


N — 

-i » va n — 


\ 2Trh? ) 


3/2 

exp 



( 21 ) 


where M = 6.7 x g is the atomic ^He mass. Be¬ 

cause of the energy-momentum conservation, the neutron 
scattering rate depends on the initial neutron state and 
on "^He atom momentum. We assume the initial neutron 
state to be given by the lowest energy level along the z- 
axis and by the momentum p|| = {px,Py\ parallel to the 
helium surface, corresponding to the total neutron energy 

K =py2m + Eo. (22) 

In typical experiments with UCN, K k, p^|/2m ^ 
10“^ eV ^ Eq 10“^^ eV. The typical initial momen¬ 
tum P of a He atom is larger than the neutron momentum 
by still more than one order of magnitude, because its av¬ 
erage kinetic energy Ebb — (3/2) kBT ^ 10“^ eV. Below 
[see Eq. (IMl) ] we will see that, as p|| —>■ 0, the neutron- 
He scattering rate rcvap remains finite, and we can neglect 
P||/P <C 1 in the calculation of this rate. 


A. Matrix elements 


The initial neutron wave function is given by Eqs. ®, 
m and ®. Since the typical neutron out-of-plane ki¬ 
netic energy after scattering by a He atom of the vapor 
is much larger than Eq, and mostly even larger than Vq 
given in Eq. 0 , the final neutron wave function is close 
to the three-dimensional plane wave with momentum p', 
normalized to one particle in the whole volume V: 

'll)' = exp (ip'r/fi) /^/V. (23) 

The initial He-atom wave function is T = exp(zPr/?i), 
and the final He wave function is 4*' = exp (fP'r/fi). The 


matrix element of the interaction potential m is given by 

Tif = J d^rV'_Lo( 2 )V'|| (i'||)V''(r) X 

J d^Rexp (r - R) 

where APtot ~ P ~ P^ + P|| — is the change of total 
momentum, R is the coordinate of the He nucleus, and 
= Cq Ai [{z — Eq/ mg) /zq] according to Eq. ([S]). 
Introducing u = zf zq and performing the integration over 
d^Ry in Eq. (1^^ using the identity 


d^R|l exp (/APtot||R||/d) = {2Tihf (APtot||), 

(25) 


one can rewrite Tjf as 

_ U(2^d)^d(2)(APtot||) 

where the remaining integral is 


(26) 


7= 1.4261 / dh Ai (m — ai) exp 


iApz u 
h/zQ 


We calculate this integral approximately by replacing the 
normalized Airy function / (ft) = 1.426lAi (m — ai) by 
a simpler form, also normalized, that is a close approx¬ 
imation, i.e. / (m) ~ exp —{u — uq)^/2 where 

uq Ri 1.5. Then, 


I = 


du f(u )exp 


/ iApz u 

V 


/■“ du 

1 

1 

"Hi 

1 

Si 

O 

to 

iAv, u 

/_oo 

2 

-1 

O 


= ^1/4^2 


exp 


iApz Up _ I f Apz y 

k/zo 2 \ h/ zq ) 


(27) 


Below we need only the square of the absolute value of 
the matrix element Tjf. The square of the d-function in 
I Tif I ^ should be treated as 


(27rd)2d(2) (APtotll) 


^( 2 ^fi) 2 d( 2 )(APtot||), (28) 


because it comes from the extra integration over the co¬ 
ordinate r^ii: f d^rqi = S. Indeed, substituting Eq. (1^ 
to the l.h.s. of Eq. (1^ we obtain 


{2Trh)'^ d^^'> (APtotll) J d^r^ii exp ^ ^APtot||U|| ^ ^ 

= {27rn)^5^^^ (APtotll) j d'uii = S i2T:hf 5^^'> (APtot||) ■ 


























5 


Substituting Eq. dUl) to Eq. (ESI) and using Eq. (l28l) . we 
obtain 


' "I V^hlz„ 

Since S/zq P-, using the identity 

1 


_ ( Apz Y 

V ) 


S (x) = lim 
£->0 


ey/TT 


exp 


we rewrite iTifl as 

(27r?i)'d(3) (APtot)/V". 


(29) 


B. Scattering rate 


The scattering rate of a neutron with initial in-plane 
momentum p|| by a He atom with initial momentum P is 
given by the square of the matrix element (l29l) integrated 
over the final momenta p' and P' of the neutron and He 
atom, respectively (Eermi’s golden rule [l^): 


27r f d^F' 

~ h J { 2 TThY 


Ed^p' 

{27rh)^ 




(30) 


Here e ~ /2M and e' = P'^/2M -I- (p — p')^ /2 to are 

the initial and final total energies of He-atom and neutron. 
The scattering rate is approximately independent of the 
initial neutron momentum p||, i.e., w-p = ri;(p||,P) ~ 
w{P), because P|| <C P can be neglected in Eq. (|5(I)) . 
We now substitute Eq. ([29|) to Eq. (1311)1 . The integration 
over p' cancels (APtot) in Eq. (1^ . where APtot ~ 
P —P' —p'. After the integration over the angle (f) between 
P and P' we obtain 


Wp = 


P'‘^dP' mU^ 

27rfi4 pp' 


(31) 


2PP' - 




Using M K. Am, from Eq. EU we obtain 


P'‘^dP' mU'^ 


Wp 


-e [8PP' - 3p2 - 5P'2] . (32) 


27rfi4 PP' 

The integrand is nonzero when the inequality 

5P'2 _ gpp/ + 3p2 ^ Q 


(33) 


is satisfied. The quadratic expression has two real roots, 
P' = P (4 ± 1) /5, so that Eq. (1331) is satisfied for 3/5 < 
P'/P < 1, thus defining the range of integration in Eq. 
((3^1) . i.e. 


P'dP'mU"^ U^Pm 8 

Xp/g 27rft4 ^ “ 27rft4 


(34) 


Einally, to obtain the total scattering rate as function of 
temperature one has to integrate Eq. (IMl) over the initial 
He-atom momentum P, weighted with the distribution 
function Np of He vapor given by Eq. Eni, 


(T) = 




-Wp 


(2Trh)^ 

P^dP^ U^8m 
{2TThf 


exp 


^ - P^/2M 


Introducing the new dimensionless variable P^/2MfcBP 
and performing the integration we obtain 


Wvap (T) 


{2MkBTf U'^Sm / /r 
{2'kK)^ 25 ^ \A:bP 


After substitution of Eq. © one obtains 


(35) 


Wvap (T) = 9.44 s 1 X (T [K])^ x exp • (36) 

Hence, rcvap (1 K) Ri 0.007 s"^ = (138 s)”^^ and the es¬ 
timated mean scattering time of a neutron in the low¬ 
est level, as determined by He vapor only, is about 138 
s at P = IK. Lowering the temperature diminishes 
the scattering rate more rapid than exponentially, e.g., 
{T = 0.8 K) Ri 7.74 x lO'^ s-^ = (21.53 min)”\ and 
Wvap {T = 0.7 K) Ri 1.6 X 10“4 s-1. The break-even with 
neutron decay is thus reached slightly above 0.8 K. 


IV. SCATTERING FROM SURFACE WAVES 


A. General information about ripplons 


A quantum of a surface wave (ripplon) with momentum 
q induces a surface deformation along the z-axis, given by 

C(i'lhi) = Cousin (qrII-Wqt) . (37) 

The dispersion relation of surface waves is given by 

Wg = — -I- q tank (qd) , (38) 

where a ~ 0.354 dyn/cm is the surface tension coefficient 
of superfluid ^He, p ~ 0.145 g/crc? is its mass density, d 
is the depth of the helium bath above a horizontal bot¬ 
tom wall, and >c^ = (g + /) p/a with an additional force 
/ oc d“4 to the van-der-Waals attraction of helium to 
the bottom wall. The ripplon amplitude in Eq. (157)) . 
normalized to one ripplon per surface area S, is given by 

[SliSlill 


For a helium bath (in fact already for a thick helium film), 
XT = \fgpja Ri 20 cm“^. The thermal ripplons with en¬ 
ergy hwq Ri k^T Ri 0.5 K have the wave number q r; 1.2 
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nm“^, for which holds q ^ x and qd ^ 1. Then the dis¬ 
persion relation of ripplons is just the dispersion of capil¬ 
lary waves: 


and 


UJg ~ 


sOg 



(40) 


(41) 


B. Interaction Hamiltonian 


To determine the influence of a periodic surface defor¬ 
mation on the neutron quantum state on the surface we 
have to separate two limits. The first, adiabatic limit ap¬ 
pears when the surface oscillates so slowly that the neu¬ 
tron wave function adjusts to the instantaneous surface 
profile. The interaction potential in this limit is found in 
Appendix A, see Eq. (EHI, and can be rewritten as 


d \ 


(42) 

where p|| = —ih ^|| and Pq — hq are the momentum 
operators of the neutron and ripplon along the surface, 
respectively. This interaction term generalizes Eq. (7) 
of Ref. [^, because it does not exclude coordinate- 
dependent surface perturbations. The expression in the 
square brackets in Eq. is just the transfer of the to¬ 
tal (neutron-fripplon) kinetic energy to the final neutron 
kinetic energy along the z-axis. 

The opposite, anti-adiabatic or diabatic limit appears 
when the surface oscillates much faster than the charac¬ 
teristic frequency of the out-of-plane neutron motion, so 
that the neutron wave function does not adjust to the in¬ 
stantaneous surface prohle. In this limit a surface wave 
affects the neutrons by creating an additional time- and 
coordinate-dependent periodic potential 


Hint — ^(i’ll , t) 


+ p<lT -Pi 

2m 


- biOo 


, , _ J Vb at 0 < z < for ^(I’lpt) > 0 

\ -Vo at C(i’lhi) < ^ < 0 for < 0 

(43) 

The ripplon amplitude, given by Eqs. (IMl) or (HD) , for 
any reasonable value of S is much less than the atomic 
scale and, even more, than the typical scale of the neutron 
wave function, given by Kq^ ~ 33 nm at z < 0 [Eq. (|5])]. 
Therefore, the potential in Eq. (1431) can be approximated 
by 


Vr{r) = I4(r||, z) Ri VbC(i’||, t)P (z) ■ (44) 


C. Crossover between adiabatic and diabatic limits 
and matrix elements 

The diabatic-adiabatic crossover, corresponding to a 
change of the ripplon-neutron interaction Hamiltonian 


from Eq. (1^ to Eq. (1421) . must take place when the rip¬ 
plon frequency ujq and the wave-vector q decrease. How¬ 
ever, the estimate of the crossover frequency tOqc and 
the description of the system in the crossover regime is 
not a trivial problem. Similar problem appears in other 
condensed-matter systems and requires a special theoret¬ 
ical study (see, e.g.. Refs, 

One may, naively, define the crossover as the region 
where the ripplon frequency becomes comparable to the 
quasi-classical bouncing frequency of a neutron in the 
ground level in z-direction, i.e. when the ratio fvjjq/Eo ~ 
1, where Eq given by Eq. 0. This corresponds to the 
ripplon frequency 

ujqc ^ Eo/h = 915 s-\ (45) 

and to the ripplon wave number 

9 c ~ = 70 cm“^ > (46) 

However, such an estimate of the diabatic-adiabatic 
crossover has an important drawback: it does not depend 
on the value Vq of the neutron potential inside helium. 
Generally, we expect that for Vq —t 0 and for non-zero q 
and LOq one can always apply Eq. (H51) . and for Vb —t 00 
one can always apply Eq. dm, which contradicts Eq. 611). 
The classical definition of the diabatic-adiabatic crossover, 
given by Eqs. (IB5I) and (IB6I) in Appendix B, has the same 
drawback. 

A rigorous analysis of the adiabatic/diabatic crossover 
should be based on the solution of the Schrodinger equa¬ 
tion for a neutron in the time-dependent potential given 
by Eqs. dSH) and (IA2|) . One may approximately deter¬ 
mine the criterion of adiabatic/diabatic crossover from 
the variational principle to minimize the neutron energy. 
This approach would be definitely correct for a time- 
independent potential. The lowest-level out-of-plane neu¬ 
tron wave function is chosen to minimize the neutron 
energy. In the adiabatic limit the energy loss is the ki¬ 
netic and gravitational energy from Eq. (1421) , while in the 
diabatic limit it is the potential energy from Eq. gal)- 
The first-order energy correction is given by the diago¬ 
nal matrix elements of these two interaction potentials, 
or more precisely by the Hamiltonian in Eq. (IA2I) . If 
these diagonal matrix elements are nonzero, their com¬ 
parison gives the crossover frequency. If these matrix el¬ 
ements vanish in the first order in one needs to cal¬ 
culate and compare the second-order corrections. Since 
JZo^^'P±oiz)di;j_o{z)/dz = 0, and / C(r||, t)d2r|| = 0, 
the first-order (in ^) diagonal matrix element of the adi¬ 
abatic Hamiltonian in Eq. (1421) vanishes. So does the 
diagonal matrix element of the diabatic Hamiltonian in 
Eq. (H51) in the first-order in if P|| ^ hq. Hence, to 
calculate the crossover frequency one needs to calculate 
the second-order energy corrections, which do not van¬ 
ish. These corrections are determined, in particular (but 
not only), by the matrix elements of the neutron-ripplon 
interaction potentials in Eqs. 611) and 61)- Therefore, 
for an estimate of the position (ripplon frequency) of the 
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adiabatic/diabatic crossover the comparison of the matrix 
elements, given below, is more accurate than just the com¬ 
parison of ripplon frequency with Eq. As we will see, the 
final result of the neutron-ripplon scattering rate is not 
sensitive to this crossover frequency, because the main 
contribution to this scattering rate comes from the rip- 
plons with energy hjjq ^ Vb ^ which corresponds 

to the far diabatic limit. 

Therefore, a rough estimate of the crossover between 
diabatic and adiabatic limits is given by the ripplon fre¬ 
quency when two interaction Hamiltonians, given by Eqs. 
(P5|) and (HU), become of the same order of magnitude. 
More precisely, we compare their matrix elements for 
the neutron transitions between lowest energy levels of 
their motion in z-direction. Thus defined, the diabatic- 
adiabatic crossover depends on Vq and meets other general 
requirements, such as the adiabatic limit for ujq,q —^ 0. 
The matrix element Tjf of the diabatic interaction poten¬ 
tial in Eq. HI for the transitions between two neutron 
states with initial wave function '4}j_ (z) and final wave 
function '0^ (z), written explicitly in Eqs. is given 

by 

^ ^oC(r||)<5(2:). 

The integral over z cancels S (z), while after substituting 
Eq. (1571) the integration over r|| gives (27r/i)^ (Aptot||) j 
where Aptot|| = fiq + p|| — py is the change of the total 
in-plane momentum of the ripplon-|-neutron system. As a 
result we obtain 

Tif = 'l^|Eo,n. (47) 

The factor 

Vii = S-^ {27rhf (48) 

is due to the in-plane part ?/;|| (r||) of the neutron wave 
function, given by Eq. and 


The first values of Qk,n are given in Table I of Ref. (2^ . 
G-g-, Qo,i = 0.09742 Qo ,2 = —0.05355 /rm“^, 

Qo ,3 = 0.03831 /rm“^, (5o,4 = —0.0304 /rm“^. Substi¬ 
tuting these values we obtain for the first levels the ratio 
bo,n/77o,n ~ 1 at 

a;,c~10^s“^ (53) 

By chance, the diabatic-adiabatic crossover condition, de¬ 
fined as Vb,ri Ho^n, is close to the value quoted in Eq. 

H 3 - 

Below, we consider mainly the neutron-ripplon interac¬ 
tion in the diabatic limit, corresponding to the ripplon 
frequency ujq > uiqc and the interaction potential given by 
Eq. HU), because it deals with much larger phase space 
of ripplons and because, as we will see later, the main 
contribution to the neutron-ripplon scattering rate comes 
from the ripplons with energy Hujq ^ Vq ^ hjjqc- 

The scattering rate rurip of a neutron with initial in¬ 
plane momentum p|| on ripplons is determined by two 
processes: the absorption and the emission of a ripplon 
with wave vector q and energy hiOq, 

^rip — a^abs 4“ '^em- (^4) 

Since for typical ^He temperatures fceT is much larger 
than the initial neutron energy K, the populations of rip¬ 
plon states with relevant energies are Nq ^ 1 or even 
Nq ^ 1. The phase volume of an absorbed ripplon is 
much larger than that of an emitted ripplon, because the 
energy of the latter is limited to the initial kinetic energy 
of the neutron, K <C A:bT. Hence, one could expect that 
the ripplon-neutron scattering rate is dominated by rip¬ 
plon absorption, so that w^ip ~ Wabs- However, because 
of a low-energy divergence of Wrip (see below), the emis¬ 
sion of low-energy ripplons with energies Hujq K may 
also be important, and we therefore consider both these 
processes. 

D. Absorption of ripplons 


Tif = 


d^i 


■^_L (z)exp iru 


Vo,n = Vo^Ogr±oi0)^±ui0) ( 49 ) 

comes from its out-of-plane part ijjj_ (z). The squared 
modulus of the matrix element in Eq. H3 follows as 


iTifl" 


{2Trhf6^^'> (Aptotii) 

S 


(50) 


where we again have used Eq. (1751) . 

The matrix elements of the adiabatic interaction poten¬ 
tial in Eq. (|42l) for ^(r||,t) given by Eq. (l37l) are 


Hk,n =V\\^Qq{Ek — En)Qk,n, (51) 

where V|| is again given by Eq. (l48l) and 

= (52) 


The absorption scattering rate Wabs of a neutron with 
initial in-plane momentum P|| in the discrete vertical level 
with energy Eq is given by Fermi’s golden rule, 


27r f NqSdPpq 

= T y {2^hY 


(27r?i)2 


^|Tif|"<5(e-e'), 


(55) 

where Pg = hq is the ripplon momentum, and pj| and 
n are the in-plane momentum and out-of-plane quantum 
number of the final neutron state, respectively. 

Nq = [exp{hwq/kBT) - 1]“^ (56) 


is the Bose distribution function of ripplons with energy 
Hujq and with zero chemical potential. The matrix element 
I Tjf I is given by Eq. (1501) and the initial total energy by 

e = Hujq +py2m + Eq. (57) 







The final energy e' = /2m+En, after using the in-plane 

momentum conservation expressed by the (5-function in 
Eq. (1501) . can be rewritten as 


p\+pI + ‘2PqP\\ coscj) 
2m 


+ En, 


(58) 


where (j) is the angle between p|| and p^. The integration 
over the component pj| of the final neutron momentum 
parallel to the surface cancels the (5-function in Eq. (l50l) . 
After substitution of Eqs. (ISO)) . (1571) and (1551) to Eq. (1551) 
we obtain 


'^abs — 


NqSpqdpqdcj) 

2T:h^ 


ElV'iO iO)r±n (0)K)eo,| 


X(5 huiq — AEn — 


pI + 2pgP|| cos (j) 
2m 


(59) 


where AE^ = — Aq ~ E^ is the change of the out-of¬ 

plane neutron energy after the ripplon absorption. The 
integration over (f) cancels the i5-function in Eq. dMl) and 
gives 


'^abs — 


NqSpqdpq |V^_L0 {0)r±n (0)M)?O,| 




■E 


a^-{b-AEr,f 


( 60 ) 

where a = p\\pq/m and b = hiOq — p^/2m.[80j We estimate 
this integral in the Appendix C. This calculation gives the 
upper estimate of Wabs [see Eqs. (ICllI) . (|C17|) and (|C21|) ] 
of 


^abs 


+ w'^ + w'^Ki 7 X 10“® X T [K] s-^ (61) 


This corresponds to a mean neutron scattering time due 
to ripplon absorption of r^jp > 4 hours even at T = 1 K. 


E. Emission of ripplons 

The rate of emission of a ripplon by a surface-state neu¬ 
tron with momentum py is given by Fermi’s golden rule, 
similar to Eq. (1551): 


_ 27r f NqSd'^Pq j" Sd'^p' ^ 

h J (2^fi)2 J (2^?l)2 E 

(62) 

Here, Pq = fiq is now the emitted-ripplon momentum, 
p' = P|| — Pg is the in-plane neutron momentum after 
emission of the ripplon, 


A^' = 1 -I- [exp (hujq/k^T) - 1] \ 


(63) 


is the ripplon population, e ~ P|/2m is the initial total 
energy, and 


, _Pl+Pq- ‘^PqP\\ cos (j) 


+ En + klOq, 


is the final total energy. The matrix element is given by 
Eq. (|5ni) . The integration over p' in Eq. (l62)) cancels the 
(5-function in the matrix element. 


SN'qPqdpq 

2TTh? 


EiV’^ ( 0 )^^:;: (o)Foeo,i 


d(j)5 


pI - ‘^PqPW cos((' 
2m 


+ En + hljJq )[64) 


The integration over the angle 4> between py and Pq in 
Eq. (IMl) is similar to that in the preceding subsection in 
Eq. dMl) and gives 


^em — 


SNqPqdpq ^ ( 0 ) ( 0 ) 




■E 


0 <j| 


— (5l -|- AEn) 


(65) 

where a = p\\Pq/m as in the previous subsection, and bi = 
Hujq + Pq/2m. The integrand is real when 0 < AEn ~ 
En < a — bi = p\\pq/m — tmjq — Pq/2m. This can be 
satisfied when a — bi >0, which for p|/2m = 100 neV 
gives Pq < pf^^ = hqf^^ with « 6.5 x 10^ cm"!. 
The maximum value of a — 6i is ^ 3 neV ^ Vq = 18.5 
neV. Hence, for the emission of ripplons, <C Vq, and we 
may use Eqs. (fT3l) . (fT51) and (IC18I) . In addition, instead of 
three intervals of parameters for the ripplon absorption, 
we only need to consider one interval. Substituting Eqs. 
(SOI), dUl) and the upper estimate of \4’'± (0)|^ < \'ipj_Q (0)|^ 
to Eq. (1651) . AV6 obtSiin 3(11 upper estimute ^ 0 ^ for 


k^Tpqdpq f dn\ip‘]_Q{0)Vo\ 


Jq 2Trhapq 

f k^Tdpq 


- (bi + En) 


( 66 ) 


^io ( 0 ) ^o| 


2—61 


yj EjidEji 


TT^h^aPq Jo 


This integral resembles the one in Eq. (IC19I) : the only 
difference is the sign of En in the denominator and, con¬ 
sequently, a different upper integration limit. We may 
give an upper estimate of this integral by replacing ^/En 
by its maximum value y/a — bi in the integrand and by 
replacing the lower limit by —bi in Eq. (1551) . This gives 


w, 


,p _ f kBTdpq |^io(0)Eo| V^^dEn 


TT^h^apq 


k^Tdpq \J 2 I 0 (0)^0I 

0 irti^apq 2gV2m 




The integral converges. Neglecting p^l2m <C EjJq and 
changing the integration variable to we finally obtain 


n"P 


fcsT |^io( 0 )Ho |^ 2 p||\3/^ 
T:b?a g\/2m 3\m/ \ a 


2m 


(67) 
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The rate of ripplon emission depends on the initial neutron 
momentum At K = p|/2to = 100 neV Eq. (1571) gives 

<P a; 2 X 10“® s"i X T [K] . (68) 

Combining Eqs. dun and (1551) we obtain an upper estimate 
for the total scattering rate of a surface neutron in the 
lowest energy level Eq by ripplons: 

<ip = '^Ihs + <m ~ 9 X 10"^ s"^ X T [K]. (69) 

This rate corresponds to a mean neutron scattering time 
due to the ripplons of Trip > 3 hours even at T = 1 K. 


V. OTHER NEUTRON SCATTERING 
PROCESSES 


A. Scattering of surface neutrons by bulk phonons 


The scattering of ultra-cold neutrons inside superfluid 
helium by bulk phonons has been studied in Ref. (d^ . 
There, two main processes were identified: (i) one-phonon 
absorption and (ii) one-phonon absorption combined with 
emission of another phonon due to the cubic term in the 
phonon Hamiltonian. The second process was found to 
dominate at low temperature, resulting in a total scatter¬ 
ing time of about Tpho = 100 s for a neutron propagating 
through liquid "^He at T = 1 K. In our case of a neu¬ 
tron above the He surface, both scattering processes are 
weakened by the factor 



V'l (0)/2 k 


1.16 X 10"'^, 


because only a small part of the neutron wave function 
penetrates into the liquid helium. Hence, for helium tem¬ 
peratures below 1 K, the neutron scattering time constant 
due to bulk phonons, Tph ~ Tpho2«:/'0l (0) > 10® s, is ex¬ 
tremely long and can safely be ignored. 


B. Scattering by surfons 

Recently, a new type of surface excitation was proposed 
( 4 ^ in addition to the ripplons, in order to explain the 
temperature dependence of the surface tension coefficient 
of liquid helium. These excitations, called surfons, are 
He atoms in a quasistationary discrete quantum energy 
level above the liquid helium surface [dsl - fs^ . The state 
is formed by the combination of the van-der-Waals at¬ 
tractive potential of the bulk helium and the hard-core 
repulsion between He atoms. Although there is so far 
only indirect experimental evidence for this type of sur¬ 
face excitation, we consider the neutron-surfon scattering 
rate to compare with the other processes. The interaction 
potential is the same as for neutron interaction with the 
helium vapor, but the surfons propagate only along the 
helium surface. Therefore, the surfon-neutron interaction 


contains an additional small factor ^ zoV'l (0) ^ 3 x 10“® 
due to a small overlap of the neutron and the surfon wave 
functions. The activation energy of the surfon has been 
obtained from fitting the temperature dependence of the 
surface tension coefficient of liquid ^He to the experimen¬ 
tal data [ 4 ^. Its value. Ago ~ A:b x 2.67 K = 3.7 x 10“^® 
erg, is significantly smaller than the evaporation energy 
—p = kB X 7.17 K of a ^He atom. Therefore, at low 
enough temperature the neutron scattering by surfons will 
exceed the scattering rate on helium vapor and must be 
considered for completeness. 

In the calculation we can neglect the initial UCN mo¬ 
mentum p as compared to the large surfon initial mo¬ 
mentum P|| ^ \/2kBT M, similarly to our treatment 
of the scattering from helium vapor in Sec. HI. We 
also assume that the surfon in-plane kinetic energy is 
not sufficient to evaporate the He atom from the sur¬ 
fon state after scattering. The vertical neutron energy 
level En may change, however, and the out-of-plane neu¬ 
tron momentum may not be conserved because the he¬ 
lium surface violates the spatial uniformity along the z- 
axis. The surfon energy consists of the excitation en¬ 
ergy Aso and of the kinetic energy Kgur — P|/2M of 
its in-plane motion. The populations of the surfon states 
are approximately given by the Boltzmann distribution. 


Nsur (P||) ~ exp 


The calcula- 


- (Ago + /2 m) /ksT 

tion is described in Appendix D and gives a very small 
upper estimate for the scattering rate of neutrons on sur¬ 
fons: 


^sur 


{MkBTf^a^Be I .2 

4 X 10-®exp ( X T3/2 [K] s-1. (70) 

V kbT ) 


kBT 


Hence, at temperatures T > 0.25 K, the neutron scatter¬ 
ing rate WgP,. by surfons is found to become much smaller 
than the scattering rate Wvap by helium vapor given in 
Eq. (1361) . However, for this and lower temperatures, the 
scattering by ripplons is dominant, <C w '^^, so that 
scattering of UCNs by surfons is negligibly small at any 
temperature. 

Thus, the total scattering rate Wtot of UCNs on the 
liquid helium surface is determined by the helium vapor 
at high temperatures T > 0.6 K, and by ripplons at low 
temperatures T < 0.6 K. It is plotted as function of tem¬ 
perature in Figs. [Hand [21 


VI. DISCUSSION AND SKETCHES OF 
EXPERIMENTAL IMPLEMENTATIONS 


The calculations presented in this paper show that, at 
temperatures below 0.7 K, the mean scattering time of 
a neutron in a gravitational quantum state above a hor¬ 
izontal flat surface of liquid helium is greater than the 
neutron beta-decay lifetime. This surface might therefore 
indeed represent an almost perfect mirror, which calls for 
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FIG. 2: The calculated total scattering rate (solid line) of 
UCNs on the liquid helium surface in inverse hours. The dot¬ 
ted and dashed lines give the contributions due to helium vapor 
and ripplons, respectively. 



T[K] 


FIG. 3: The logarithm with base 10 of the calculated total 
scattering rate of UCNs on the liquid helium surface. 

experimental demonstration and applications. The sys¬ 
tem could offer excellent possibilities not only to study 
the quantum states represented in Fig. [T]but also serve as 
a sensitive probe for detection of tiny energy transfers due 
to helium-intrinsic or external perturbations. The appli¬ 
cation having motivated this theoretical work is a high- 
precision study of the level scheme of the neutron in the 
gravitational potential above the liquid mirror, giving ac¬ 
cess to short-range, gravitation-like interactions between 
the neutron and the mirror. A motivation from an exper¬ 
imental point of view has been current work on new UCN 
sources at the ILL [H, which involves cooling many 
liters of ultrapure, superffuid helium below 0.7 K. The de¬ 
velopment had been started at the TU Munich 
and builds on theoretical work by Golub and Pendlebury 
on superthermal UCN production via down-scattering of 
cold neutrons in superfluid helium 


The scattering rate of neutrons on helium vapor de¬ 
creases stronger than exponentially while lowering the 
temperature [oc exp (—const/T), see Eq. (|55)) ]. Al¬ 

ready for 0.8 K it is calculated to be smaller than the 
neutron decay rate. An upper estimate of the scattering 
rate due to ripplons at 0.8 K is found to be by one or¬ 
der of magnitude lower than the rate due to the vapor. 
Owing to its linear temperature dependence it will be¬ 
come the dominant contribution below about 0.6 K [see 
Eq. (I69p ]. however at a level already 50 times below the 
neutron decay rate. While all processes calculated here 
should be insignificant for precision studies of the level 
scheme, experiments involving storage of neutrons with 
energies up to the cutoff set by the neutron optical poten¬ 
tial barrier Vq of the superfluid helium [see Eq. ©] may 
have different requirements. In this respect it seems help¬ 
ful that the main contribution to the neutron scattering 
rate at such low temperatures is due to the low-energy 
part of the ripplon spectrum and thus will dominantly 
lead only to transitions between nearby neutron quantum 
states. Energy transfers of a few peV are however usu¬ 
ally insufficient to cause a neutron to penetrate through 
the liquid helium and thus leave the system. Therefore, 
at 0.6 K the mean escape time of a UCN with initial ki¬ 
netic energy K < Vq could be longer than the neutron 
beta-decay lifetime by several orders of magnitude. This 
makes an experimental set-up using a liquid helium sur¬ 
face a strong candidate for a nearly loss-free neutron con¬ 
tainer and has indeed been proposed to be applied in a 
neutron lifetime experiment [^, [^. Eor highest relia¬ 
bility, measurements should nonetheless be performed at 
different temperatures and the container preferentially be 
filled with a neutron spectrum with a gap between its up¬ 
per cut-off and Vq. It should also be noted that the value 
Vq = 18.5 neV for superfluid helium is small compared 
to > 100 neV for conventional materials used for neu¬ 
tron bottles. Counting statistics might therefore become 
a limiting issue. Still, a neutron lifetime measurement 
employing a trap involving a horizontal surface of super¬ 
fluid helium seems an interesting complement to projects 
employing magnetic neutron traps. While these possess 
typical trapping potentials for low-field-seeking neutrons 
in the range (50 — 120) neV and completely avoid any wall 
collisions of truly trapped neutrons, other systematic ef¬ 
fects such as marginally trapped neutrons and depolariza¬ 
tion need to be carefully addressed [60l - l^ (see also Refs. 
[H, 1^ for recent reviews and discussion of the neutron 
lifetime problem). 

Turning to the question how to populate and detect the 
neutron quantum states above a superfiuid-helium mirror 
(called ’Take” in the sequel) one first notes that, in con¬ 
trast to a solid mirror as employed in previous and ongoing 
experiments, one has to confine the liquid and deal with 
the presence of a meniscus at the border of the helium con¬ 
tainer (the ’’coast” of the lake). In the traditional ’’flow- 
through” scheme of current experiments neutrons enter a 
mirror table with an absorbing ceiling from one side and 
are detected on the other side. This might also be tech- 









11 


nically feasible for the liquid mirror, where neutrons will 
then however have to enter and exit the lake through a 
thin, weakly absorbing foil with low (or negative) neutron 
optical potential. Application of magnetic fields might of¬ 
fer more attractive, novel experimental possibilities which 
we sketch below. 

Considering first a technique for state population, we 
note that one may employ a magnetic field gradient for 
deceleration of neutrons moving from above towards the 
horizontal mirror located at z = 0. Neutrons with mag¬ 
netic moment in a magnetic field with modulus B have 
a potential energy of ± B, with sign depending on the 
spin state with respect to the field direction. The upper 
(lower) sign refers to those neutrons which become re¬ 
pelled (attracted) by a positive gradient of magnetic field 
modulus. They are correspondingly called low-field (high- 
held) seekers. Note recent experimental demonstrations of 
trapped high-held seeking UCNs A magnetic held 

modulus B (z) = Cz, for instance with constant (posi¬ 
tive) gradient C > mg/ Ri 1.66 T/m overcompensates 
gravitation for the high-held seeking neutrons. Those with 
vertical kinetic energy E±^ = {\g,^\C — mg) h at height h 
will thus have lost this energy entirely when arriving at 
the mirror. This situation is analog to a neutron rising in 
the earth’s gravitational held to its apogee (after which it 
will fall back down). If alternatively, one wants to dec- 
celerate the low-held seeking neutrons, the gradient has 
to be inverted and hence the strongest held needs to be 
located at the surface. While the neutron is close to the 
lowest point of its trajectory, the magnetic held gradient 
nearby the mirror needs to be switched off. For limitation 
of the spatial region where the held needs to be provided, 
one may employ a vertical, straight neutron guide, ending 
above and close to the mirror. A circular absorber with 
a central hole for the neutron-feeding pipe and mounted 
with variable distance of some tens of /rm above the lake 
may serve for preparation of low neutron quantum states 
as used in the hrst experiments [T^. An obvious beneht 
of the magnetic population method is a possible neutron 
detection acceptance angle of the full 2tt. Using a helium 
lake offers the additional advantage that the presumably 
nearly perfect mirror can be made much larger than the 
quartz mirrors employed in current experiments. 

Compared with a how-through experiment the sensitiv¬ 
ity of the energy state determination may be drastically 
improved using lateral UCN trapping prior to detection, 
ideally for many hundreds of seconds. This might be possi¬ 
ble using a magnetic fence, consisting of a multipolar mag¬ 
netic held arrangement similar to the system described in 
Ref. [zH . For our purpose the multipole of high order has 
to be oriented such as to provide only held components 
within the plane dehned by the helium surface, i.e., with 
current-carrying rods parallel to the gravity held. This 
will keep low-held seeking neutrons away from the liquid 
meniscus at the container wall, thus keeping the verti¬ 
cal neutron state unaffected. Note that, if one populates 
the lake with high-held seeking neutrons, their spin has 
to be hipped after arrival at the surface and prior to ar¬ 


rival in the region of strong multipolar held, which can be 
done, e.g., using standard magnetic resonance techniques. 
Detection can still be done through the side walls of the 
helium container, requiring switching-oh of the magnetic 
helds (or better a spin hip to turn the trapped low-held 
seeking neutrons into high-held seekers to accelerate them 
through gaps in the magnetic fence). Alternatively, one 
may let them rise back to the entrance of the neutron 
guide by switching on again the magnetic held gradient 
used for lake population. 

Next, we discuss a possible further improvement of 
the lake population technique, noting hrst that a verti¬ 
cal straight and specularly rehecting guide with constant 
cross section does not mix the vertical with the horizontal 
components of neutron velocity. As a result, the closer 
a typical neutron approaches the mirror, the larger will 
be the number of rehections per length unit of the guide. 
Even a small non-specularity in the rehection will then be¬ 
come an issue. In addition, a typical neutron from a typi¬ 
cal UCN source will, after removal of its kinetic energy in 
the vertical direction, still have a typical hnal speed par¬ 
allel to the surface of several meters per second. A lower 
neutron speed would however be benehcial for both the 
how-through mode and for trapping. For the former it 
increases the time a neutron spends on the mirror, while 
for the latter a lower magnetic held strength is sufficient 
for lateral UCN conhnement on the lake. 

A non-imaging neutron optical device proposed by 
Hickerson and Filippone ohers an interesting remedy 
for the aforementioned dehciencies of a straight guide 
[t^. They describe a compound parabolic concentrator 
(CPC) for neutrons rising from a Lambertian horizontal 
disk source upwards against the gravitational held. Its 
neutron-optical properties are based on the ” neutron foun¬ 
tain” valid for constant force helds along the symme¬ 
try axis of a parabolic rehecting surface. Using a constant 
magnetic held gradient that overcompensates the ehect 
of gravitation, neutrons approaching the lake from above 
will experience a constant deceleration a = C/m — g. We 
may thus apply the CPC inverted in space with the neu¬ 
tron source (an aperture with radius R) located at height 
h above the lake. According to the formulas given in Ref. 
[t^ , neutrons starting there at time t = 0 and with speed 
vq will, after a time T = vr/ a (where vr = \/vq + 2aR) 
and with typically fewer than two rehections, arrive within 
a narrow band of heights 0 < z < R above the horizon¬ 
tal mirror. The spread of total kinetic energy within the 
ensemble of UCN is then only AE r; maR/3 and inde¬ 
pendent of uq. After switching off the held gradient, the 
neutrons close to the mirror will thus move with much re¬ 
duced lateral velocities compared to the traditional beam 
method. Hence, even without lateral trapping, state pop¬ 
ulation via a CPC will lead to much increased interaction 
times with the mirror and a corresponding gain in accu¬ 
racy of measurements. A CPC will be most benehcial at 
a pulsed UCN source, preferably in combination with a 
rebunching technique as demonstrated in Ref. [t^, and 
hence works best in combination with a UCN trapping 
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experiment. We note that the very low lateral neutron 
velocities allow for quite modest magnetic trapping fields, 
which makes it easy to provide large openings for neutron 
detection in the magnetic fence. Obviously, a CPC with 
a magnetic deceleration system might also be used for a 
sufficiently large conventional mirror. 

In summary, this paper has given positive answers 
concerning necessary prerequisites for application of a 
superfluid-helium mirror for study and application of neu¬ 
tron quantum levels in the Earth gravity field. Further 
investigations will be needed to address questions on neu¬ 
tron manipulation, e.g., if transitions between levels can 
be induced by vibrations of the helium surface in a con¬ 
trolled way. It might also be worthwhile to consider fur¬ 
ther possibilities to create a flat mirror of large surface, 
such as ’’Fomblin” oil (a fluorinated, organic compound 
with low neutron absorption, already tested as part of 
an optical system for a new neutron charge measurement 
0) or liquid or solid neon. 

P. G. thanks A. M. Dyugaev for useful discussions. The 
work was supported by RFBR grant #13-02-00178. 


Appendix A: Derivation of the nentron-ripplon 
interaction in the adiabatic approximation 


The Schrodinger equation for a neutron is given by 

H'ljj {r,t) = ihdip {r,t)/dt, (Al) 

where the Hamiltonian 

H = K + V = — -I- mgz + VqO \—z + ^(r||, t)] (A2) 

contains the neutron kinetic energy K = — and 

the potential energy V = mgz + VqO [—z + CCru, 0] ■ The 
latter contains the effects of the Earth’s gravitational field 
and the potential wall due to the liquid helium, as shown 
in Fig. [T] The difference from Eq. (|3]) is that the liq¬ 
uid helium has now a time- and space-periodic boundary 
^(r|| ,t) given by Eq. (l37l) . The difference between Eq. 
(|A2|) and Eq. (4) from Ref. [1^ is that the surface has 
now a periodic spatial dependence. 

The adiabatic adjustment of the neutron wave func¬ 
tion to the new surface profile means that the neu¬ 
tron wave function, in first approximation, adiabatically 
shifts in z-direction by the length ^(r||,t): #(r,t) —>• 

# (r -I- ^(r||, t)£, t), where z is the unitary vector in z- 
direction. This shift can be written via the translation 
(z-shift) operator 

Tz (C) = exp [^(r||,t)a/az] = exp [i^{r\\,t)p,/h] . 

Its action on the wave function is 

{z) = tl’^{z + C)- 


We also define a new wave function 


# (r -f t) = # (r, t) = {£) # (r, t ), 

which after substitution into Eq. m gives a new 
Schrodinger equation for 'll; (r, t): 

Ht (0 ^ (r, t) = ihd (f, (0 # (r, t)) /dt. (A3) 

The action of the shift operator on the potential energy 
function V (r) is given by 

r (r) t (0 = t io V{r- e(r||, t)z) , (A4) 


while for the commutator with kinetic-energy operator 
K = —T?Kl2m we have 


KT^ (0 -f^{0K = -— (\;/2e^irvt)d/dz _ ^a’rn,t)d/dz^2 

'm V 


2p||Pg + pI 

2m 




-PtY-pfi 


2m 




(A5) 


where py = —y and Pq = hq are the neutron and 
the ripplon momentum operators along the surface, re¬ 
spectively. The time-dependence of ^(ry , t) also gives an 
additional term on the r.h.s. of Eq. (IA3l) : 


“dt 
= ih 


= ih 

f^d^p (r, t) 




d 


dt 


+ hxjq^{r\\,t)—T^ip{r,t). (A6) 


Combining Eqs. (IA3I) and (IMI) we obtain a new 
Schrodinger equation. 


+ Hint - ihd/dt'^pj {r,t) = 0, (A7) 


where Hq is given by Eq. o and the interaction term is 
given by 


Hint — C(j'|| #) 


2p||P9 + pI 
2m 



d 


(A8) 


Appendix B: Crossover between adiabatic and 
diabatic limits in classical physics 

For a classical particle above the surface in the limit 
Vo —> oo the crossover between diabatic and adiabatic 
limits occurs when the maximal acceleration of the helium 
surface d^£,q/dP = due to its oscillatory motion, 

becomes equal to the free fall acceleration g, 

= 9- (Bl) 

The classical amplitude of the surface oscillations with 
wave vector q differs from in Eq. (ITO by the square 
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root of the Bose distribution function Nq given by Eq. 

(ElDiizl 

~ ioq\JkBT/hljJq. (B2) 

In addition, in Eq. (ITO depends on the surface S, 
which must be defined. In the formulas for the neutron 
scattering rate by ripplons this surface-dependence is un¬ 
physical and does not occur explicitely, because the S- 
dependence of the ripplon amplitude in Eq. dun is com¬ 
pensated by the 5'-dependence of the ripplon density of 
states (see below). Similarly, the total mean square am¬ 
plitude of thermal surface oscillations at any point r|| is 
given by the sum over all q-vectors, 

h)> = E^q= / 


and the surface area S drops out. More generally, if we are 
interested in the surface waves with the wave number q in 
some interval Aq^Agy, then we sum all ripplon modes in 
the phase volume SAq^Aqy, and the surface area S again 
drops out from the total (C^)- In the estimate (IBII) for 
the diabatic-adiabatic crossover the surface S is defined 
by the area Sn of the neutron wave function along the 
surface, which corresponds to the momentum smearing 
Aq^Aqy S'“h 

For a lower estimate of leading to an upper estimate 
of the ripplon scattering rate, we take the minimal possible 
S given by the square of the wave length: S'min ~ {2n/qf. 
Then, substituting it to Eq. dUD, we have 


^max 


hq^ \ ^ / hl2 \ 

2-k V\/^/ 


(B3) 


For q = 1 p,m ^ = 10^ cm ^ this formula gives ~ 7 x 
10“® nm, which is much less than q~^. The corresponding 



(B4) 


is also much less than and we can apply the usual 
surface wave description. Substituting S’min ~ [27:/q) to 
Eqs. (HD), (IBIII and (155)) gives 


a 


= w3/2 
^qc 


( kBTq^^ 

1^2 ( 2^)2 


which, using Eq. (go)), gives the lowest estimate for the 
crossover frequency 


[^qc — 



6.6 X 10’^ s-i 

{T[K\f^ 


(B5) 


This frequency corresponds to the neutron energy (at T = 
1 K) hujqc ~ 7 X 10“^^ erg Ri 0.44 neV > Eq and to 


qc = 0.56 fim ^ > >c. (B6) 


Hence, in the diabatic limit q > qc, and we can always use 
the ripplon dispersion given by Eq. (HOI). 

Another condition of the classical adiabatic limit is that 
the curvature of the surface = 9 ^^™“ is less than the 
curvature of the neutron trajectory due to the parabolic 
free-fall motion where r:|| is the neutron 

velocity along the surface. Taking a UCN kinetic energy 
of ATii = 100 neV, corresponding to Uy = 2Ar||/TO = 19 
m^/s^, we can check that the condition < g/v‘^\ 

is fulfilled at q = qc = 0.56 Hence, the condition 

ujq < ujqc given by Eq. (IBs!) ensures the classical adiabatic 
limit. 


Appendix C: Calculations for the neutron scattering 
rate due to ripplons 

In this section we evaluate the integral in Eq. dSD or 
O, which gives the neutron scattering rate by ripplons. 
The integration over pq and n in Eqs. (l59l) can be sep¬ 
arated into several regions, given by different limits of 
the ratio Pq/p\\ and of the difference AEn — Vq. For 
AEn < Vo the final neutron vertical state belongs to a dis¬ 
crete energy spectrum, approximately given by Eq. (HD. 
For AEn > Vo the final neutron vertical state belongs to 
the continuous energy spectrum and can approximately 
be taken as a plane wave. 

In the region pq ^ py the initial neutron kinetic en¬ 
ergy is negligible and, for the majority of the scattering 
events, the change of the neutron out-of-plane kinetic en¬ 
ergy AEn > Vo- The integral in Eq. dSD is evaluated in 
this limit in Appendix Cl below. 

In the region of small ripplon momentum, pq < py, 
studied in Appendix C2, the angle (j) between the ini¬ 
tial neutron and ripplon momenta is important for the 
out-of-plane energy transfer AEn, and the scattering rate 
depends on the initial neutron momentum py. Depend¬ 
ing on the sign of the difference AEn — Vo, this region 
is split into two. For AEn < Vo the final neutron state 
belongs to the discrete spectrum and is described by the 
formulas in Sec. H. For AEn > Vo the final vertical neu¬ 
tron state belongs to the continuous spectrum and can be 
approximated by Eqs. (IC2p and (IC3I) . 


1. Absorption of thermal (high-energy) ripplons 

In this subsection we consider the region of large mo¬ 
menta Pq ^ p|| = \/2Km contributing to the integral in 
Eq. (I59p . Let us assume the in-plane kinetic energy K of 
ultra-cold neutrons being less than K„ = 100 neV, which 
corresponds to a maximal initial neutron wave number 
a* = 7*11*/^ = 7 X 10^ cm“^ = 0.07 nm“^ and to a max¬ 
imal neutron velocity riy* = \/2Kjm = 4.4 m/s. For 
q = q* the ripplon energy, according to Eq. (go]), is given 
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by 

fiwq* = hujq {hq = p||*) ~ 600 neV ~ 7 mK Vb, 

and the ripplon velocity is Vq^ = 3wg»/2g* = 20 m/s [TSj. 
If a ripplon with such a high energy is absorbed, the fi¬ 
nal out-of-plane neutron energy ^ hwq is much higher 
than the potential barrier Vq = 18.5 neV. It is then reason¬ 
able to take the final out-of-plane neutron wave function 
as a plane wave, 

i’±n {z) ~ exp {ip'^zlh) (C2) 

Accordingly, the neutron out-of-plane energy can be ap¬ 
proximated by the free-particle quadratic dispersion 

En^p^l2m, (C3) 


dlol) and jMl) to Eq. jCTj, and introducing the dimen¬ 
sionless variable ^ = fnOq/k^T = h^Ja/pq^/'^/k^T, for 

which q = (^(kBT^p/a/h^ , we obtain 

- («/<„„)■/’ ' 
(C 8 ) 

where Cmin = fuZq*/kBT is given by Eq. (jCTI) and = 
= (2m)^ (a/hp)^ /^bE 1. The integration 

in Eq. (IC8I) diverges as Cm[n^ ^t the lower limit, and the 
main part of the integral comes from this divergence: 


IV’uo (0) bbl^ ^/mkBT 
3TrH‘^aV2y^C~ 


(C9) 


where p/ is the component of the final neutron momentum 
perpendicular to the surface. The sum over out-of-plane 
neutron wave number n in Eq. (1591) then becomes an in¬ 
tegral over p'^: 

^ / Pnip'z)dp'^, (C4) 

n 

where the one-dimensional neutron density of states is 
given by [t^ 

Pn (P'z) ~ L^/2TTh. (C5) 

Eor scattering by thermal ripplons with q > q^, the 
initial neutron energy K « p|/2m and the momentum 
P|| < pq can be neglected. Eqs. (1571) and (1551) then sim¬ 
plify to 

e Ri hwq] e' K, (p^ +p'z) 1‘^m. (C6) 

Using Eq. (1551) . we rewrite Eq. dMD as (the lower index 
”>” means large ripplon energy): 


Substituting the cutoff = tvjjq^/kBT given by Eq. 
(EU and other numerical values to Eq. (IC9I) . we obtain 
the contribution to the neutron scattering rate from the 
high-energy ripplons with hq > p||*: 

!(;> R! 1.7 X 10”® s-i X r[K]. (CIO) 

At smaller ripplon energy, i.e. at /iq < p|p, the integral 
in Eq. (l59l) must be estimated without the approximation 
in Eqs. (IC2p - ()C6l) . In the next subsection we show that 
Eq. ()C9p overestimates the integral in Eq. (l59l) for hq < 
P||, especially for hujq < Vb where the infrared divergence 
disappears. 

At AT, —0, when the cutoff given by g* = ?i/p||* is too 
small, the infrared divergence in Eq. (IC9|) must be cut off 
Cmin ~ kb/^BT, because the approximation given by 
Eqs. (IC2I) - (IC5I) is not valid for lower ripplon energies, for 
which the neutron state after the absorption still belongs 
to the discrete spectrum along the z-axis. A rough esti¬ 
mate of the absorption rate of high-energy ripplons with 
kijJq > Vb can be obtained for small initial neutron ener¬ 
gies AT < Vb by substituting Cmin ~ Vb/feeT to Eq. ()C9I) : 


w> 




(0) K)| NqPqdpq 

2'Kh? 




IV'uo (O)^ol^ mNq^dq 

2Trhy/pa ^2mhujq — h^q"^ 


(C7) 


The square root in the denominator is real at 2mhijq = 
2mh\foiJp(^l'^ > h'^q^, which gives q < 4qo = 

[2m/K) ajp ~ 2.5 nm“^ and corresponds to the rip¬ 
plon energy hiOq < hujqmax ~ h{2m/h) (a/p) = 2 x 

10-16 gj.g = 1.25 X 10-4 eV R IcB X 1.5 K > feT. Above 
this energy the simple absorption of a ripplon by a UCN 
in a surface state is impossible because of the conserva¬ 
tion laws for energy and momentum. Substituting Eqs. 


R w> (Cmin ~ Eo/^bT) r IQ-"^ s'^ X T [K] . (Cll) 

This estimate gives a neutron mean scattering time 
i/w'^ R 27 hours, which is much greater than the in¬ 
trinsic neutron lifetime. 


2. Upper estimate of the absorption rate of 
low-energy ripplons 

For /cbT ^ hujq the ripplon population is given by Nq r 
kBT/hujq. Substituting Eqs. (15(71) and (HTl) to Eq. (1551) we 
obtain 


w< = 


kBTdpg ^ \tp±Q ( 0 ) tp±n ( 0 ) 7 /)| 

Jo 4 ^a 2 -(&-An)^ 


(C12) 
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a. Transitions to a continuous neutron spectrum 

In this subsection we consider the case of final neutron 
energies En > Vo above the potential barrier and thus 
belonging to a continuous spectrum. We may then apply 
the approximation given by Eqs. (IC2I) - (IC5I) and rewrite 
Eq. (ICT^ as 


r*"" fcsTdPg |V^_Lo (0) Vbl^ 

47rtep2 2TTh 

where the integral 


(C13) 


1= f ^ r+'’ Vml2E„dEn 

\J a? — {h — /2mf' \J a? — {b — En)^ 

(CM) 

For b > a we may give an upper estimate of this integral: 

T T — dEn TT^/m 

^ --vmLu 7:rf^-75%- 

(C15) 

The corresponding upper estimate of Eq. (jClSp is 


_up |'i/^_i_Q (0) Vq I I Pmax 

XV ^ ^ /_ m 


4:TTh?'a^2Vo 


Pmin 


(C16) 


The interval of integration Vq < En < 6 + a in Eq. (IC14I) 
is nonzero for 6 + a ~ fwjq + p\\Pq/m > Vq, which for 
p|/2m = 100 neV corresponds to g > qmin = Pmin/d ~ 
4 X 10"^ cm“^. Substituting also '0 _lo (0) ~ 0.236 cm“^/^ 
and Pniax = P\\* to Eq. (IC16E we obtain 


w<P Ri 7.3 X 10"® In {q^/qmin) x T [K] r; 2 x lO"® x T [K]. 

(C17) 


for En < Vq. Using Eq. (IC18|) . we rewrite Eq. (IC12I) for 
AEn < Vq as 


/■P— kBTdpq \^lo (O)Uof y^dEn 

^ Jo nfi^apq ngV^ Jo 

(C19) 

This integral converges, with main contributions from 
En ^ Vq. The upper estimate of this integral can 
be obtained by replacing y/E^ by in the integrand 
and by extending the integration region from (0, Vq) to 
{a — b, a + b). This gives an integral over En similar to 
Eq. lIM : 



Since Vo > En and the integrand in Eq. (IC19I) is real for 
b—a < En < 6+a, the region of integration over En in Eq. 
(IC19I) is nonzero if 6 — a r; hiOq — p\\Pq/m < Vq, which for 
p|/2to = 100 neV corresponds to Pq < Pmax ~ fi^maxyo 
with qmaxvo = 1-57 X 10® cm“^. On the other hand, 
AEn < b + a can reach Vq ii b + a Hujq + p\\Pq/m > 
Vq. For Py/2TO = 100 neV this gives Pq > Pminyo ~ 

fi^minyo with gininyo = 3.8 x 10"^ cm“b For pq < Pminyo 
the logarithmic divergence disappears. Hence, using Eq. 
(IC20|) we obtain an upper estimate of w^: 

R! 4 X 10“® X T [K] s-C (C21) 


b. Transitions to the discrete neutron levels 


Appendix D: Scattering of neutrons by surfons 


In this subsection we consider the case of final neu¬ 
tron energies in the interval 0 < AEn ^ Vo below the 
potential barrier and approximately given by Eqs. (USD 
and USD. Since Vo ^ Eo, the sum over n in Eq. (IC12I) 
still includes many terms and can be approximated by an 
integration over n for n ^ 1. Eqs. (US and ([13 give 
En ~ mgzo (37rn/2)^^^, which can be rewritten as 


Stt \ mgzo J 


2{AEnf^ 

Sirgh 



and gives 


dn 

dEn 



(C18) 


We also use that |'0j_„ (0)| < tj}j_Q (0), and for an upper 
estimate of Eq. (IC12I) we replace |V'ln(9)| by IV'.lo (0)1 


1. Matrix element 

A surfon, being a ^He atom on the surface energy level 
[ 4 ^ [ 4 ^ , interacts with a neutron via the potential given 
in Eq. O- The matrix elements of neutron-He interaction 
in Eqs. (123,® assumes that the He wave function is a 
plane wave also along the 2 ;-axis, which is not the case for 
the surfons. Therefore, in this subsection we derive the 
neutron-surfon matrix element in the way similar to that 
in Sec. HIA. 

The surfon wave function is given by a product 'I' = 
(z) 4>|| (ry). The parallel-to-surface surfon wave func¬ 
tion is a plane wave: (ry) = exp (zP|| r||/?i) . The 

perpendicular-to-surface surfon wave function 'I'x (z) was 
analyzed in Refs. [49l - [5i| and shown to be localized above 
the surface on a height ^ 0.4 nm <C zo [8l[. Hence, for our 
calculation we may take (z) r: S (z). The initial and fi¬ 
nal surfon wave functions differ only by the initial and final 
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momentum, Py and Py, respectively. On the other hand, 
the neutron out-of-plane wave function strongly changes 
due to the scattering on a surfon, and for the majority of 
events it gets transferred from the discrete to the contin¬ 
uous spectrum. Hence, the final neutron wave function is 
given by Eq. (1251) . The matrix element of the interaction 
potential (ED is given by 


Introducing new dimensionless integration variables Xp = 
P|P/P| and Up = Pj_/P||, we rewrite this as 


U‘^mP\\ r |'*/'io(0)| dxp dtjp 
^ ^ \J 64xp - (3 -I- 5xp + 


(D6) 


Ti{ 


U 


WV'ii (i'll) V'' (r) X 


This may be further transformed to 


{z) exp 


*(P||-P')R|| 


<5^^^ (r - R) 


where APtot ~ P — P' — p' is the change of total momen¬ 
tum and R is the surfon coordinate. Below we need only 
the square of the absolute value of the matrix element Tif. 
The square of the d-function in |Tif|^ should be treated 
using Eq. (1251) . Then instead of Eq. (122)) we obtain 


iTif = C/^V'io (0) (APtotii) IV. (D2) 


Wp 


U'^mP\\ 

U^mP\\ 


Iv^io (0)| dxp dyp 
\/(Vp - dl) {b2 - y^) 
IV’lo {0)\dxp dyp 

~ ®i) (“2 - Xp) ’ 


where the solutions of square-root equations are 
3 ^ SyjXp 5xp 


(D7) 

(D8) 


(D9) 


2. Scattering rate 


The scattering rate of a neutron by a surfon with initial 
momentum Py is given by the square of the matrix ele¬ 
ment (1221) integrated over the final momenta p' and Py of 
the neutron and surfon, respectively (Eermi’s golden rule 


^ f f Vd^p' 

h J {2'KhY J {2TTKf‘ 


iTifl'd (£-£')■ (D3) 


and 


17T8,/l-20y2-20y2 

ai.2 =--. (DIO) 

The integral over Xp in Eq. (IDSp gives tt for any 02 > ai, 
however, the integrand is real only for some values of y^. 
The maximum value of &2 in Eq. (1 d 21) is = 0.2 at 
Xp = 16/25. Since 0 < for an upper esimate of 

the integral in Eq. ()D6p we may take 


Here e ~ P|/2M and e' = P'^^/2M+p''^/2m are the initial 
and final total energies of surfon -|- neutron, respectively. 
We now substitute them and Eq. (ID2I) to Eq. (ID3I) . The 
integration over py cancels (APtot||) in the matrix 
element in Eq. (ID2I) . where APtot|| ~ P|| ~ Py ~ 


i,“P 


C/^mPy 

P^mPy 


f |'*A±o ((^)| dXp 2 



- ai )(02 
2Try/0l2 

5 ■ 


Xp) 


(Dll) 


U^lPlo ( 0)1 

[ ^\\^^\\ dPx [ 


1 2tt 2TTh J 


[ d(j)5 is — s'). 
-n J 


(D4) 


The total scattering rate is given by the integral over all 
initial surfon states with corresponding populations: 


Substituting e and e', and integrating over the angle (j) 
between P y and P y, we obtain 



Now we use M = Am to simplify this expression: 


,Ji„(0)| J ” ‘^( k^T )• 


Substituting Eq. (|Dlip and performing the integration we 
obtain 


w 


up 

sur 


V'io (0) 


P^m\/l.67r /—Aso\ 




Wp = 


{U'^m\ip\Q[0 )\/A-k’^P^) P||dP||dp(j 


p/2p2 _ 


P|f/8 - P|f/8 - (P|f + P|f + p'l) /2]' 


Substituting U from Eq. ED we obtain the scattering rate 
given by Eq. (EQD, which is negligibly small. 
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